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TECHNICAL MEMORANDUM X- 64963 


A TETHER TENSION CONTROL LAW FOR TETHERED SUBSATELLITES 
DEPLOYED ALONG LOCAL VERTICAL 


SECTION I. INTRODUCTION 


A. Historical Development of Tetherline 

The use of tethers has been proposed In the past for a multitude of pur- 
poses. One significant study effort was directed towards the retrieval of 
stranded astronauts by using tetherlines [1, 2, 3], The approach was to "throw 
a buoy" from a rescue vehicle towards an astronaut. The buoy would be 
attached to the astronaut and then reeled in. 

Reference 1 contains a derivation of relative motion equations. Also 
included are nomogr£q)hs which enable one to determine the Initial ejection 
velocity relative to the main spacecraft and the constant tether tension required 
to cause the trajectoiy to pass through the desired point. The scheme, as 
presented, is open loop; that is, there is no control over perturbations from 
the desired trajectory. 

The significant contribution of reference 2 was a descxiption of the 
problems associated with retrieval, especially die terminal phase, when 
relative velocities and angular rates can become quite large. Here again, the 
control was open loop, thus small perturbations in the system at the beginning 
of the retrieval can cause large deviations at the terminal phase. 

Reference 3 presented a concept for a reel mechanism for controlling 
line tension and a "fish pole" to aid in preventing tanking. Rather stringent 
tolerances were required to assure that the control system could control the 
trajectory as desired. A closed loop scheme was implied, but the control law 
was not discussed. 

Another body of work was directed towards station keeping requirements. 
Reference 4 describes a scheme designed to tether a modifiled lunar module and 
Apollo telescope mount (LM-ATM) to an orbital workshop (OWS). This scheme 


required periodic tugs on Ihe tetherline to cause the LM-ATM to move about 
the OWS. A repetitive* somewhat egg shaped trajectory relative to the OWS 
could be obtained by accurately timing and metering tether tension during the 
tugging operation. The scheme presented safety hazards and required constant 
attention by the crew and was judged to be cm unatti'active flight configuration. 

Manned flight e:q>erience with tethers was obtained on Gemini XI and 
Xn. Reference 5 describes these experiences. The Gemini XI flight demon- 
strated a rotating configuration with the Gemini spacecraft tied to an Agena by 
100 ft of dacron webbing and spinning at 0. 9 deg/s. The Gemini Xn fli^t 
demonstrated gravity gradient ctqpture. Insertion into the stable equilibrimn 
was difficult because of thruster failures, but peak amplitudes of less than 60 
deg from local vertical were observed. 

The Radio Astronomy Explorer (RAE) flight e:q)erience is of interest [6] . 
The RAE uses a combination of magnetic torquing, long (750 ft) booms, and 
hysteresis damping to bring the veMcle close to stea<fy state equilibrium. To 
achieve rapid damping, the booms were deployed at opportune times in the 
spacecraft attitude oscillation. The effect was to decrease the attitude rate by 
Increasing the inertia of the spacecraft. Subsequent deployment of a hysteresis 
damper nq>idly reduced libration amplitudes to stea^ state. Although this 
e}q>erience does not involve tethers, the rather long booms exhibit some of the 
same characteristics of tetheilines and some of the same techniques can be 
used to damp tetherline motion. 


B. Objectives of This Study 

All of the previously mentioned projects, with the exception of the RAE 
and Gemini Xn flight, involve the use of tetherlines (or long booms) in dynamic 
modes of operation. The objective of this study, however, is to design a con- 
trol law which will cause the motions (swingir ' and stretching) of a tetherline 
to be damped to the steady state equilibrium. It 'will be shown that structurcl 
damping in the tetherline will cause motion of a tethered subsatellite to be 
damped. Normally, this damping is neaiiy inconsequential, but an active 
system will be shown which enhances the damping to the point where even large 
amplitude oscillations in the orbital plane can bt damped to nearly zero in two 
cydes. 


The equations of motion presented are planar to simplify the presenta- 
tion, but the contid law has been studied in a three degree-of-freedom relative 
motion orbital simulation. The effects of atmospheric drag have been lumped 
at the subsatellite and the tetherline is assumed to be massless. 
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SECTION II. EQUATIONS OF MOTION 


Figure 1 portrays a small subsatellite tethered to die main satellite, in 
this case the orbiter. It is assumed that the orbiter mass is much larger than 
that of the subsatellite and the orbit is circular so the orbit rate, ^ , is con- 
stant. From reference 1 the equations of motion for the subsatellite rdative to 
the main body, considering motion only in the orbital plane and a massless 
tetherline, are: 


f(0 + 0)2 + 


(la) 
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- A”2) sin 0 

e 

(lb) 
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where 

c = tedierline coefficient of viscous damping 

K = tetherline spring constant 

i - unstretched tetherline length 
o 

H = operator whose value is 1 if f a f ^ , 0 if t < 

ri = radius of the orbiter orbit 
T = tether tension . 
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Equation (la) is i'eferred to as the length equation and (lb) as the swing 
equation. Inspection of equation ( lb) shows that 0 = 0 , tt are equilibrium solu- 
tions; that is, 0 = 0 when ) = 0 = 0 and 0 = 0 or it. At the equilibrium solu- 
tion: 

9 = 0 . = 1 - afj + - 10(7} ■ • • (2«) 

9 = ir , A"* = 1 + 3^^ ♦ 6(^J + w(^) • • • (2b) 

Settii^ 0 = 0 and taking the first 2 terms for in equation ( la) 

gives: 
for 0 = 0, 


= t(i)‘ * -i^) + <j ( 3 .) 


for 0 = 7T , 


}• = l(i)‘ . i«r,[2i . 3(ijj . H[|(i - V * . (3b) 


Dropping liie j term and assuming positive tension reduces equation 
( 3) to a linear differential equation: 


.. C • /k •A K 

-M,' • 


( 4 ) 


The parameters for damping ratio, ^ , and natural frequency, , for 
the stretch equation are: 
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Note that the usual natoral frequency of a spring mass system is modified 
by a gravity gradient term, 3^^ . This inq>lies that Ihe spring constant term, 
K/M 2 , must be larger than the gravity term if steady state equilibrium is to 
exist. The final value of t is given by 


t 

ss 


Ma 


f 

o 



( 6 ) 


t 


I 

i 

1 


i 

! 


Equation (6) n^lects higher order terms due to gravity gradient. An error of 
approximatdy 2% exists in the stretch due to this omission which is acceptable 
for most analyses. 

The natural frequency and damping ratio of the swing equation, ( lb) , 
can be approximated in a similar fashion assuming constant i . For e dose 
to 0 or 7T, equation (lb) is approximated by 


•• X • *0 i • 

fl = - 2j 0 - 30^0 - 2 j 0 


(?) 


The damping ratio and natural frequencj^ are 

(jj ss sfz 0 ( 8a) 

n 



n 


(8b) 


Inspection shows that the swinging frequency is independent of £ and 
that the swinging motion is damped during deployment where i > 0 , The damp- 
ing can be significant as can be seen from an illustration. Consider a deploy- 
ment rate of 5 m/s, a length of 10 km, and ^ = 0.00 1185 rad/s. The damping 
ratio is 0. 24. 

Maintenance of a fixed length or retraction of a tetherline causes the 
damping term to disappear or change polarity (indicating instability). Addi- 
tional means of control must be employed to enable satisfactory performance 
for the latter two operations. 


SECTION 1 1 1. DEVELOPMENT OF A TETHER CONTROL LAW 


Quite often, a system which has coupling from one mode of vibration to 
another can be stabilized by adding damping to one mode only. The effective- 
ness of this procedure depends, among other things, on the d^ree of coupling 
between the modes. In the case of the tethered subsatellite, viscous damping 
in the tetherline aids in damping the swinging motion. The viscous damping of 
materials which could be used for tetherlines is very low* However, a tether 
reel located in the orbiter can be controlled in a manner which simulates any 
damping or spring constant. Implementation of such a control law requires 
sensing tether tension, rate of line deployment, and length of line deployed. 

A computer calculates the tension as a function of dre rate and line length 
required to simulate a spring with the desired characteristics. The computer 
then commands the reel motor to produce the required tension. The tension at 
the subsatellite end of the tetherline approximates the tension at the orbiter end 
if the tether mass is ne^igible. If not n^ligible, the mass of the line causes 
delay and change of amplitude of tension which can be predicted and compensated 
by the computer. 

Computer simulation results will now be shown which illustrate the 
dEunping of initial conditions for natural damping obtained from tetherline 
materials and controlled damping obtained from a reel mechanism control law. 
Consider a tetherline system with the following characteristics: 


0 = 0.001185 rad/s 
K = 0.2 N/m 
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Mj = 100 kg 

C = 0.1 N-sVm 

f = 4 km 
o 


A trajectory of this system in a rotating local vertical coordinate frame 
is shown in Figure 2. The subsatellite is initia' :ed at 4 km viiich is approxi- 
mately tile unstretched tetiierline length. The tetiier tension was zero whenever 
the length was less than approximately 4005 m as indicated by the dotted trajec- 
tory. Length increases positively in a downward direction. A small amount of 
damp in g is evident due mainly to structural damping which is 1. 1% of critical 
damping. Figure 3 shows a plot of the swing angle, 6 , versus time and n^- 
ligifale damping of the swinging motion can be observed. 

From the preceding example, additional control is required if initial 
conditions are to be damped. It is evident that littie coiqpling exists from the 
swinging motion to the stretching motion. This is due in part to the wide 
difference in the resonant frequencies of the two vibrations as seen by comparing 
equation (5a) to equation (8a). Equation (8a) shows that the swing frequency is 
a function of only the orbital rate and no adjustment is possible, but equation 
(5a) includes the spring constant, K, and mass of the subsatellite vhich can be 
adjusted. For a given subsatellite mass, the tetherline spring constant can be 
selected to yield a wide range of stretch natural frequencies. The tetherline 
material can be selected to give the desired spring constant or it is much pre- 
ferred to design a reel mechanism control law to give an effective spring con- 
stant as required. The effective spring constant programmed into the reel 
mechanism controller is in series with that of tiie tetherline, and the tetiierline 
spring constant can be made negligible compared to that of the reel mechanism. 
Since it is desired to make the stretch frequency equal to the swing frequency, 
the required spring constant can be foimd by equating equations (6a) and (8a) 
and solving for Kj: 

Ki = 60“ M2 . (9) 


Similarly, an effective coefficient of viscous damping can be progranuned 
into the reel mechanism. Equation ( 5b) can be used to determine the coefficient 
required for a given damping ratio: 



( 10 ) 



The reel mechanisin control law used has essentially the same form as 
equation (id): 


T 


K,i + c,i + K„i 
1 1 2 c 


( 11 ) 


Instead of » \\hicn is the initial unstretched length of the tetherline* 
i denotes the commanded lengdi. The factor K. is found from an equation 

C 6 

similar to equation (4): 



( 12 ) 


The steady state solution of ( 12) is 



Forcing ^ = gives the value required for K^; 

K2 = Mj - 30*) . (14) 


The control law* equation (11) , commands tethexiine tension as a func- 
tion of amount of line played out, rate of line deployment or retrieval and 
commanded length. Thxis, three sensors are required to implement the control 
law. 


The effectiveness of the control law can be demonstrated using the follow- 
ing characteristics: 


8 



M2 = 100 kg 

i = 4100 m 
c 

i = initial value of i = 4000 m 
o 

? = 1.0 

0 = 0.001185 rad/s . 


Using equation (9) , K| = 8. 43 E-4 N/m; using equation ( 10) . Cj = 0. 410 
N-s/m; and using equation ( 14) , K 2 = 4. 22 E-4 N/m. A tetherline system with 
this control law has a trajectory as shown in Figure 4. The swing angle as a 
function of time is shown in Figure 5. This demonstrates the greatly improved 
results obtainable by controlling the tetherline. 


SECTION IV. DEPLOYMENT OF TETWERLINES 


One of the applications of tetherlines required deployment of tetherlines 
100 km long [7] . The quickest way to show the feasibility of deploying such 
long tetherlines is to use the control law developed in Section m. The sub- 
satellite could be given an initial ejection velocity and the commanded length in 
the control law could be adjusted as required to obtain '*n acceptable trajectoiy. 

A problem arose which caused the subsatellite to free-fly. If the com- 
manded length is too large compared to the present length and velocity, the 
tension commanded fay the control law is negative which is an impossible (aera- 
tion to perform on a flexible tetherline. An immediate fix was applied which 
consisted of incrementing the commanded length. Table 1 gives the commanded 
length and increments to achieve a deployment of 100 km. The resulting 
trajectory is shown in Figure 6. Plots of tetlier acceleration and deployment 
rate are shown in Figures 7 and 8, respectively. 

One can envision other control laws uhlch could effect the dcplo 3 rment 
in a more optimal fashion, eliminating the jerky d^loyment rates vhen the 
commanded length is changed. This job, left for future investigations, is beyond 
the scope of this preliminary feasibility study. 
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TABLE 1. DEPLOYMENT SEQUENCE 



• Aerodynamic drag is neglected 

• i n^lects drag and h^her order gravity gradient terms 

c 

• Initial conditions; i = 1 m, * =0.5 m/s, 0 = 180*, 0=0 

o o o o 


• Subsatellite is lowered in steps by changing i at times indicated 

c 


Time (hr) 

0 

0.74 

1.47 

2. 18 

2. 90 

3.61 

4.32 

5.03 

(km) 

0.42 

0.74 

1.29 

2.26 

3. 95 

6.92 

12.1 

21.2 

Time (hr) 

5.75 

6.46 

7.78 

7. 89 

8.60 

9.31 

10.02 

10.73 

(km) 

30 

40 

50 


70 

80 

90 

100 


Note: The times given are the times at which 0 was observed to 
pass through 180 deg in a positive direction. 


SECTION V. RETRIEVAL OF TEMRLINES 


The concluding remarks of Section n pointed out the negative damping of 
the swinging motion vhich results when tetherlines are retrieved. This problem 
can be overcome by at least two approaches. One possible solution is to accom- 
plish the retrieval in steps, letting the control law damp the motion between the 
steps. Another scheme is to retrieve the subsatellite slow enough so that the 
positive damping in the control law is still sufficient to overcome the negative 
damping in the swing equation. The latter scheme might be more sensitive to 
initial conditions H.nn the former, but should be less time consmning. 

An example of the first approach is shown in Figure 9. The trajectory 
shown is the result of commanding the control law to retrieve the tetherline to 
50 km in length from an initial 100 km. Damping of the swinging motion induced 
is accomplished in approximately 1-1/2 hours. 
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An example of the second scheme was not performed. Again, the study 
was limited to a determination of feasibility, so only one promising approach 
was pursued. 


SECTION VI. INCLUSION OF AERODYNAMICS 


The very long tetherline Eqpplication proposed by reference 7 involved 
operation at altitudes of 90-120 km vdiex'e aerodynamic drag can not be ne^ected. 
Accordingly, the equations in Section n were modified to include die effects of 
aerodynamic drag on the subsatellitc and the teiherline itself. It is felt that a 
good approximation is to lump the total drag of the entire tetherline system at 
the subsatellite. 

Figure 10 depicts the aeroctynamic model of the tetherline sys^mi. If 
the tetherline is along local vertical so that 0 = ir , the altitude of the sub- 
satellite, A , is 
s 

A = A - f , (15) 


where A^ is the altitude of the orbiter and f is tetherline length. The drig 
on an incremental length of tetherline is given by 

AFd = ”P V^CjjDAf , (16) 

where 

p = atmospheric density at the altitude of the segment 
V = velocity of the tetherline segment relative to the atmosphere 
Cjj = drag coefficient of the tetherline 

D = diameter of the tetherline 


Ai = length of the segment of tetherline. 
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Figure 11 presents a plot of the drag exerted on each kilometer of a 
topical tetherline along local veiUcal using data from the U. S. Standard 
Atmosphere* 1962. The plot shows that the lower few kilometers of tetherline 
is )^ere the majority of the drag is concentrated. This phenomenon justifies 
the assumption that the total system drag can be liunped at the subsatellite. 

The density* p * varies with the altitude ( or length) so the total drag on 
the tetherline must be found by integrating equation ( 16) : 



(17) 


Since the tetherline is generally not along local vertical (d # ir)* the 

altitude of the subsatellite* A * can be approximated given i * 0 * and A : 

s o 


A « A - f cos (0 - v) 
so 


(18) 


Another consequence of 0 # is that the cross sectional area of the tetherline 
presented to the atmospheric drag is reduced by cos (0 - tt). Equation (17) 
becomes 

i=A -icos(0-7r) 

Fd * I Cd D / ° p cos (0 - ff) dl . (19) 


A very good exponential curve fit can be found for the drag force v^ich 
sliminates Ihe need to evaluate the integral in equation ( 19) . Assuming the 
orbiter is at 200 km and A^ is given in km* 

f=A icos(0-7r) 

f ° pcos (0 - jt) df « 10** EXP (18.057 - 0.45602 A + 0.001495 A J) . 

i=0 ® ® 

( 20 ) 
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The above curve fit is suitable for subsatellite altitudes ran g in g from 150 km 
dovm to 80 km. For hi^er altitudes» the aerodynamic drag is negligible so the 
drag force is set to zero. 

The acceleration of llie subsatellite caused by the drag force is given by 


^d “ M + 


( 21 ) 


where M is the subsatellite mass and M. is the mass of the tetherline. The 

S X 

total drag acceleration can be divided into two components: 




( 22 ) 


T , = - X cos 0 
a a 


(23) 


where i , is directed along the tetherline and T , is perpendicular to the 
d d 

tetherline. The angular acceleration with respect to the rotating frame caused 

byT^ is 


0 = — 
®d I 


(24) 


These accelerations caused by aerodynamic drag can now be added to equations 
(la) and (lo): 


i = i(0 + hy + 0* rj (1 - A”*) cos 0 - ^ ^ ^ ^d 


0* = - 2 -^ (^ + 0) " ^ (1 “ 0 + *0 d , 


(25b) 
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An example has been simulated using the following parameters: 


A 

= 200 km 

<P = 

0.001185 

o 




: 7304 m/s 

r< = 

6571 km 

V = 


c.. 

= 2.2 

Ki = 

1. 651 E-3 

D 


K2 = 

9. 436 E-4 

M 

= 168 kg 



s 


c, = 

3. 98 

“i 

= 24 kg 

D = 

3.66 X 10"® 


The initial conditions for this case are: 

i = 0 , i = 70 km , f =80 km 

c 

9 = 0 f 6 = ir . 

At intervals of 2650 s, the value of i was increased by 10 km to a final value 

c 

of 110 km. The trajectory of this simulation is shown in Figure 12. 


SECTION VI I. CONCLUSIONS AND REC0M^".ENDATI0NS 


The deployment of long tetherlines (100 km) has been shown to be 
promising from a mechanics viewpoint, but several mechanics and <fynamics 
questions remain op^. The out~c. plane motion was not analyzed in this study, 
and only a few simulations of the dc yment of tetherlines were made on a 
three degree-of-freedom relative motion simulation (8) . The results of these 
few simulations d emons trated that out-of-plane motion was bounded during 
deployment and steacfy state operation with more out-of-plane motion occurring 
in sl mi ilationa of polar orbits than for 28. 5 deg inclined orbits. The out-of- 
plane motion for inclined orbits was excited by out-of-plane components of the 
relative wind caused by the rotating atmosphere. 

Humpin g of inplane motion excited during tetherline retrieval was simu- 
lated, but it is e 3 q>ected that damping of out-of-plane motion will be difflcult. 
The analysis should be extended to include out-of-plane motion so that suitable 
tedhniques can be studied for control. 
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Additionally, a modal analysis should be performed on a distributed mass 
model of the tetherline. This work is required to verify that tetbeiiine d^mamics 
is not a ^/roblem. 

These comments are relative to the feasibility of tethered subsatellites 
from a mechanics viewpoint Additional work also needs to be done in the areas 
of 1hermo(fynamics, materials, orbiter interaction, communication, tracking, 
subsatellite design, ejqperiment acconunodation, etc. It is hoped that the 
analysis and design of the control law presented in this report are a contribution 
to the tethered subsatellite technology. 



Figure 1. Small subsatellite tethered to the main satellite. 
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DIRECTION OF FLIGHT 



Figure 2. Trajectory of tethered subsatellite with no active control, 










Figure 5. Swing angle, 9, as a function of time for a tether system with active control. 




Figure 6. Deployment Trajectory, 






Figure 7. Telher accderation. 
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Figure 10. Aerodynamic modd of tetherline system. 
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Figure 11. Aerodynamic drag on each kilometer of tefherline versus altitude 
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Figure 12. Trajectory witb atmospheric drag. 
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